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Abstract 

The paper establishes weak convergence in (7(0,1] of normalized stochastic processes, generated by 
Toeplitz type quadratic functionals of a continuous time Gaussian stationary process, exhibiting long- 
range dependence. Both central and non-central functional limit theorems are obtained. 
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1 Introduction 

Let {X{t), t S R} be a centered real-valued stationary Gaussian process with spectral density f{x) and 
covariance function r(t), that is, r{t) = f{t) = /(a;) dx, t G R. We are interested in describing the 

limit (as T —> oo) of the following process, generated by Toeplitz type quadratic functionals of the process 
X{t): 

pTt pTt 

Qrit) = / / g{u — v)X{u)X{v) dudv, tG[0,l], (1-1) 

Jo Jo 

where 

g[t) = f e^^*g{x)dx, tGR, (1.2) 

Jr 

is the Fourier transform of some integrable even function g{x), a: G R. We will refer to g{x) and to its Fourier 
transform g{t) as a generating function and generating kernel for the process Qrit), respectively. 

The limit of the process dm is completely determined by the spectral density /(x) (or covariance function 
r{t)) and the generating function g{x) (or generating kernel g{t)), and depending on their properties, the 
limit can be either Gaussian (that is, Qrit) with an appropriate normalization obeys a central limit theorem), 
or non-Gaussian. The following two questions arise naturally: 

(a) Under what conditions on f{x) (resp. r{t)) and g{x) (resp. 'g{t)) will the limit be Gaussian? 

(b) Describe the limit process, if it is non-Gaussian. 

Similar questions were considered by Fox and Taqqu [5], Ginovyan and Sahakyan [5], and Terrin and 
Taqqu [25] in the discrete time case. 

Here we work in continuous time, and establish weak convergence in G[0,1] of the process (11.11) . The 
limit processes can be Gaussian or non-Gaussian. The limit non-Gaussian process is identical to that of in 
the discrete time case, obtained in [25] . 

But first some brief history. The question (a) goes back to the classical monograph by Grenander and 
Szego m, where the problem was considered for discrete time processes, as an application of the authors’ 
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theory of the asymptotic behavior of the trace of products of truncated Toeplitz matrices (see [TB], p. 
217-219). Later the question (a) was studied by Ibragimov [17] and Rosenblatt [23], in connection to the 
statistical estimation of the spectral function F{x) and covariance function r(t), respectively. Since 1986, 
there has been a renewed interest in both questions (a) and (b), related to the statistical inferences for long 
memory processes (see, e.g., Avram [T], Fox and Taqqu [B], Ginovyan and Sahakyan [5], Ginovyan et al. m, 
Giraitis et al. m, Giraitis and Surgailis m, Giraitis and Taqqu Terrin and Taqqu |26j , Taniguchi and 
Kakizawa [24| . and references therein). In particular, Avram [1], Fox and Taqqu [ 6 |, Giraitis and Surgailis 
[SI, Ginovyan and Sahakyan [^ have obtained sufficient conditions for the Toeplitz type quadratic forms 
(5 t( 1) to obey the central limit theorem (GLT), when the model X(t) is a discrete time process. 

For continuous time processes the question (a) was studied in Ibragimov m ( in connection to the 
statistical estimation of the spectral function), Ginovyan [ 7 ] [5] , Ginovyan and Sahakyan m and Ginovyan 
et al. where sufficient conditions in terms of f{x) and g{x) ensuring central limit theorems for quadratic 
functionals ( 5 t( 1 ) have been obtained. 

The rest of the paper is organized as follows. In Section [5] we state the main results of this paper 
iTheorems 12.11 - 12.4|) . In Section [3] we prove a number of preliminary lemmas that are used in the proofs of 
the main results. Section |T| contains the proofs of the main results. 

Throughout the paper the letters C and c with or without indices will denote positive constants whose 
values can change from line to line. 


2 The Main Results 


In this section we state our main results. Throughout the paper we assume that /, 5 G L^(]R), and with no 
loss of generality, that g > 0 (see [TU], [T3]l. 

We first examine the case of central limit theorems, and consider the following standard normalized 
version of CU: 

QT{t)-.= T-^/^{QT{t)-nQTm, te[0,l]. (2.1) 

Our first result , which is an extension of Theorem 1 of [10] , involves the convergence of finite-dimensional 
distributions of the process Qxit) to that of a standard Brownian motion. 


Theorem 2.1. Assume that the spectral density f{x) and the generating function g{x) satisfy the following 
conditions: 


f-gGL\R)nL^ 


and 


nOO 

1E[Q^(1)] —>■ IGtt^ / f'^{x)g^{x)dx as T ^ 00 . 

J 00 

Then we have the following convergenee of finite-dimensional distributions 


( 2 . 2 ) 

(2.3) 


Qrit) crB{t), 

where Qxit) is as in h2.1]) . B(t) is a standard Brownian motion, and 

/ OO 

P{x)g^{x)dx. (2.4) 

-00 

To extend the convergence of finite-dimensional distributions in Theorem l2.1l to the weak convergence in 
the space Cp, 1], we impose an additional condition on the underlying Gaussian process X(t) and on the 
generating function g. It is convenient to impose this condition in the time domain, that is, on the covariance 
function r := f and the generating kernel a := g. The following condition is an analog of the assumption in 
Theorem 2.3 of [15] : 

r(-) G L^’(R), a(-) G L'^(K) for some p,q>l, i(2.5) 

p q 2 
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Remark 2.1. In fact under (12.21) . the condition (12.51) is sufficient for the convergence in (12.31) . Indeed, let 
P = p/{p — 1) be the Holder conjugate of p and let q = q/{q — 1) be the Holder conjugate of q. Since 
1 < < 2, one has by the Hausdorff-Young inequality and (12.51) that ||/||p < Cp||r||p, H^Hg < Cg||a||q, and 

hence 

i + l = 2-i-i<l/2. 

p q p q 

Then the convergence in (|2.3I) follows from the proof of Theorem 3 from m- Note that a similar assertion 
in the discrete time case was established in |13| . 

Remark 2.2. Observe that condition (12.51) is fulfilled if the functions r(t) and a{t) satisfy the following: 
there exist constants C > 0, a* and /3*, such that 

|r(t)|<C(lA|tr*-i), \a{t)\<C{lA\tf-^), (2.6) 

where 0 < a*,P* < 1/2 and a* + /3* < 1/2. Indeed, to see this, note first that r(-), a(-) G L°°(R). Then 
one can choose p,q > I such that p{a* — 1) < —1 and g(/3* — 1) < —1, which entails that r(-) G LP(R) and 
a(-) G L'^(R). Since 1/p + 1/q < 2 — a* — /3* and 2 — a* — /3* > 3/2, one can further choose p,q to satisfy 
1/p + l/q > 3/2. 

The next results, two functional central limit theorems, extend Theorems 1 and 5 of m to weak conver¬ 
gence in the space (^[0,1] of the stochastic process Qr/t) to a standard Brownian motion. 

Theorem 2.2. Let the spectral density f{x) and the generating funetion g{x) satisfy condition 12.‘A) . Let the 
covariance function r(t) and the generating kernel a{t) satisfy condition 12. 5\) . Then we have the following 
weak convergence in Cp, 1]; 

Qrit) => crB{t), 

where Qrit) is as in 12. 1\) . a is as in /I2.4D, and B{t) is a standard Brownian motion. 


Recall that a function u{x), x G K, is called slowly varying at 0 if it is non-negative and for any t > 0 


lim 


i{xf) 
0 u(x) 


1 . 


Let 5'Vb(]R) be the class of slowly varying at zero functions u(x), a: G M, satisfying the following conditions: 
for some a > 0, u(x) is bounded on [—a, a], lima;_>o w(a:) = 0, u(x) = u(—x) and 0 < u(x) < u(y) for 
0 < X < y < a. An example of a function belonging to S'Vb(R) is u(x) = |ln with 7 > 0 and a = 1. 

Theorem 2.3. Assume that the functions f and g are integrable on M and bounded outside any neighborhood 
of the origin, and satisfy for some a > 0 

fix) < |a:|"“Li(a:), |5(a:)| < \x\~<^L 2 ix), x G [-a, a] (2.7) 


for some a < 1, (3 < 1 with a + ft < 1/2, where Liix) and ^2(2^) are slowly varying at zero functions 
satisfying 

L, G SVoiR), x-^^^+^^L.ix) G L^[-a, a], i = 1,2. (2.8) 

Let, in addition, the covariance function r(t) and the generating kernel a{t) satisfy condition 12.5\) . Then 
we have the following weak convergence in C[0, 1]; 

Qrit) => crB{t), 

where Qrit) is as in 12.1\} . a is as in and B{t) is a standard Brownian motion. 

Remark 2.3. The conditions a < 1 and /3 < 1 ensure that the Fourier transforms of / and g are well 
defined. Observe that when a > 0 the process {X{t),t G Z} may exhibit long-range dependence. We also 
allow here a + ft to assume the critical value 1/2. 
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Remark 2.4. The assumptions f ■ g € T^(R), f,g € L°°(R\ [—a, a]) and (12.81) imply that f ■ g € T^(R), so 
that in (12.41) is finite. 

Remark 2.5. One may wonder, why, in Theorem 12.31 we suppose that Li{x) and L 2 {x) belong to S'Vb(R.) 
instead of merely being slowly varying at zero. This is done in order to deal with the critical case a+j3 = 1/2. 
Suppose that we are away from this critical case, namely, f{x) = |a;|““/i(x) and g{x) = |a;|“^/2(a;), where 
a + /3 < 1/2, and li(x) and l 2 {x) are slowly varying at zero functions. Assume also that f{x) and g{x) are 
integrable and bounded on (— oo, —a) U (a, +oo) for any a > 0. We claim that Theorem 12.31 applies. Indeed, 
choose a' > a, P' > P with a' + P' < 1/2. Write f{x) = \x\^li{x), where 6 = a' — a > 0. Since 

li{x) is slowly varying, when \x\ is small enough, for some e £ (0,5) we have Then one 

can bound \x\^~^ by c|ln|x||~^ £ S'Vo(R) for small |x| < 1. Hence one has when |a::| < 1 is small enough, 
f{x) < \x\~‘^' |ln . Similarly, when |x| < 1 is small enough, one has g{x) < \x\~^' |ln |a::| . All 

the assumptions in Theorem 12.31 are now readily checked with a, P replaced by a' and /3', respectively. 

Now we state a non-central limit theorem in the continuous time case. Let the spectral density / and the 
generating function g satisfy 

f{x) = \x\~°‘Li{x) and g{x) = \x\~^L 2 {x), a; £ R, a < I, P < 1, (2-9) 

with slowly varying at zero functions Li(a;) and L 2 {x) such that J-^\x\~°‘Li{x)dx < oo and |a:| -^L2 {x)dx < 
00. We assume in addition that the functions Li{x) and L 2 {x) satisfy the following condition, called Potter’s 
bound (see formula (2.3.5)): for any e > 0 there exists a constant C > 0 so that if T is large enough, 
then 

+ i = l,2. (2.10) 

Note that a sufficient condition for (12.101) to hold is that Li{x) and L 2 {x) are bounded on intervals [a,oo) 
for any a > 0, which is the case for the slowly varying functions in Theorem 12.31 

Now we are interested in the limit process of the following normalized version of the process Qxit) given 
by (11.11) . with / and g as in (12.91) : 


ZT{t) := 


1 

T°‘+f^Li{l/T)L2{l/T) 


{QT{t)-nQT{m- 


( 2 . 11 ) 


Theorem 2.4. Let / and g be as in with a < 1, P < 1 and slowly varying at zero functions Li{x) and 
L 2 {x) satisfying h2.10\) . and let Zxit) be as in \2.11\) . Then for a P > 1/2, we have the following weak 
convergence in the space Cp, 1].' 

Zrit) Z(t), 

where the limit process Z (t) is given by 


Z{t)= [ Ht{xi,X2)Widxi)Widx2), (2.12) 


with 


Htixi,X2) = \xiX2\ 




' git{x2 — u) _ 2 ^' 

i{xi + u) 


i{x 2 — u) 


(2.13) 


where W{-) is a complex Gaussian random measure with Lebesgue control measure, and the double prime in 
the integral i2.12\) indicates that the integration excludes the diagonals xi = 2 iX 2 - 


Remark 2.6. Comparing Theorem l2.4l and Theorem 1 of [5^, we see that the limit process Z{t) is the same 
both for continuous and discrete time models. 

Remark 2.7. Denoting by Pt and P the measures generated in Cp, 1] by the processes Zxit) and Z{t) 
given by (12.111) and (12.121) . respectively. Theorem 12.41 can be restated as follows: under the conditions of 
Theorem 12.41 the measure Pt converges weakly in Cp, 1] to the measure P as T — >■ oo. A similar assertion 
can be stated for Theorems 12.21 and 12.31 


4 






















It is worth noting that although the statement of our Theorem l2.4l is similar to that of Theorem 1 of |25) . 
the proof is different and simpler, and does not use the hard analysis of [25], although some technical results 
of [IS] are stated in lemmas and used in the proofs. Our approach in the CLT case fTheorems 12.11 - I2.3L 
uses the method developed in m, which itself is based on an approximation of the trace of the product of 
truncated Toeplitz operators. For the non-CLT case fTheorem 12.4p . we use the integral representation of 
the underlying process and properties of Wiener-Ito integrals. 


3 Preliminaries 

In this section we state a number of lemmas which will be used in the proof of the theorems. The following 
result extends Lemma 9 of m- 

Lemma 3.1. Let Y(t) be a centered stationary Gaussian process with spectral density fvix) S L^(R)nL^(R). 
Consider the normalized process: 


Lrit) := ^ (^£\^iu)du-M Y^u)du ^ . 


Then we have the following convergence of finite-dimensional distributions: 

/ OO 

fY{x)dx, 

-OO 


(3.1) 


(3.2) 


where B{t) is standard Brownian motion. 

Remark 3.1. Observe that the normalized processes Qrit) and Lxft), given by (12.111 and (13.1|) . can be 
expressed by double Wiener-Ito integrals (see, e.g., the proof of Lemma 13.91 belowl. In our proofs we will 
use the following fact about weak convergence of multiple Wiener-Ito integrals: given the convergence of the 
covariance, the multivariate convergence to a Gaussian vector is implied by the univariate convergence of 
each component (see [21], Proposition 2). 

Proof of Lem,m,a \S. 1[ For a fixed t, the univariate convergence in distribution 


LtIJ:) N^OCay) as T —oo 


follows from Lemma 9 of [10]. To show (1^ . in view of Remark 13.11 and Proposition 2 of [2T] , it remains 
to show that the covariance structure of Lyit) converges to that of (JYB{t). Specifically, it suffices to show 
that for any 0 < s < t, 


IE 


(Lrit) — Lt{s)Y 


ay ■ {t — s) 


as 


OO. 


(3.3) 


Indeed, using the fact that for a Gaussian vector (Gi,G2) we have Gov(Gi,G2) = 2[Gov(Gi, G2)]^, and 
letting ry{u) = e^^^ fY{x)dx be the covariance function of Y{t), we can write 


IE 


{Lrit) — Lt{s)) =2{t — s) 


CT{t-s 


1 - 


-T(t-s) \ ~ s) 


r{u)du. 


Since fvix) G L^(K.), the Fourier transform ry(u) G L^(R) as well. So by the Dominated Gonvergence 
Theorem and Parseval-Plancherel’s identity, we have as T —>■ oo 


IE 


(LT(t) — Lt(s)) 2(t — s) / ry(u)du = 47r(t — s) / fy(x)dx = ay(t — s). 


(3.4) 


□ 
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We now discuss some results which allow one to reduce the general quadratic functional in Theorem 12.11 
to a special quadratic functional introduced in Lemma r3.il 

By Theorem 16.7.2 from [18) . the underlying process X{t) admits a moving average representation; 

/ CO pCO 

a{t — s)B{ds) with / \d{t)\‘^dt < oo, (3.5) 

-OO j —CO 


where -B(t) is a standard Brownian motion, and d(t) is such that its inverse Fourier transform a(x) satisfies 
f(x) = 27r|a(a:)p. Assuming the conditions (12.21) and (12.31) . we set 

b(x) = (27^)l/2o(x)(5(a;))^/^ 


and observe that the function b(x) is then in L^(R) due to condition (|2.2I) . Consider the stationary process 


V(t) = 



b(t — s)B(ds) 


(3.6) 


constructed using the Fourier transform b(t) of b(x) and the same Brownian motion B{t) as in (13.51) . The 
process Y{t) has spectral density (see [TU], equation (4.7)) 

fvix) = 2Trf{x)g{x). (3.7) 

We have the following approximation result which immediately follows from Lemma 10 of [lOj . 

Lemma 3.2. Let Qrit) be as in i2.1\) and let Lrit) be as in hS.l]) with Y(t) constructed as in it 5*. 61) . Then 
under the conditions i2.2\) and Ii2.3\) . for any t > 0, we have 

lim Var[(57’(t) — Lrit)] = 0. 

T-»-oo 

The following lemma is a straightforward adaptation of Lemma 4.2 of m for functions defined on K. 
Lemma 3.3. If pj > 1, j = 1,... ,k, where k > 2 and 'p ~ ^ ~ then 


P k 

/ \fi{xi)... fk-i{xk-i)fk{xi + ... +Xk-i)\dxi ...dxk <Y\ WfjWpi- 

j=i 


The following lemma will be used to establish tightness in the space C'[0,1] in Theorem 12.21 

Lemma 3.4. Let the covariance function r(t) and the generating kernel a(t) satisfy condition h2.5\) . and let 
QT{t) be as in i2.1\) . Then for all 0 < s < t < 1 and T > 0, there exists a constant C > 0, such that 


IE 


IQT(t) — Qr(s)P <C{t — s). 


(3.8) 


Proof. For convenience we use the Wick product notation: : X{u)X{v) := X{u)X{v) — IE [A('u)A(u)]. So 


for 0 < s < t < 1, we can write 


rTt pTt 


/■Ts rTs 


Qrit) — Qt{s) = 




a{u — v) : X{u)X{v) : dudv — 


10 JO 
rTt rTt 


to JO 

rTs rTt 


a(u — v) : X{u)X(v) : dudv j 
dudv 


^ Pi I Pi I (2 Pi s Pi I 

= —j= / / a{u — v) : X{u)X{v) : dudv-\ — - 1 = / / a{u — v) : X{u)X{v) : 

vT Jts Jts VT Jo Jts 

: = A(s, t, T) + B{s, t, T). 
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Now we estimate B{s, t, T) (the function A(s, t, T) can be estimated similarly). We have by Theorem 3.9 of 
[19] that 

^ pT s I'T s pTt 

Wi\B‘^{s,t^T)\ = — / dui / dvi / du 2 / dv 2 a{ui — Vi)a{u 2 —V 2 )^{'-X{ui)X(vi) w X{u 2 )X{v 2 )‘■) 
^ Jo Jts Jo Jts 

^ pT s pTt pT s pTt 

= — / dui / dvi / du 2 / dv 2 a{ui - vi)a{u 2 - V 2 )[r{ui - U 2 )r{vi - V 2 ) + r{ui - V 2 )r{vi - U 2 )] 

^ Jo Jts Jo Jts 

: = Bi{s,t,T) + B 2 {s,t,T). 


By the change of variables Xi = ui — ui, X 2 = V 2 — U 2 , X 3 = U 2 — ui, X 4 = V 2 , and noting that r(-) and o(-) 
are even functions, we have 


4 r 

Bi{s,t,T) < — / dx 4 / \a{xi)a{x 2 )r{x 3 )r{xi + X 2 + X 3 )\dxidx 2 dx 3 . 
d Jts dR3 

Since |r(f)| < r(0), we have r(') £ L°°(K.). We also have r(-) £ LP(K) by condition (12.5L where l/p+ 1/g > 
3/2. The L^’-interpolation theorem states that if a function is in and with 0 < pi < p 2 < 00 , then it 
is in LP , Pi < p' < p 2 . By the L^’-interpolation theorem, one can choose p' > p such that r(-) £ (M) and 


1 1 1 1 o , . 1 1 

—7“t“ —— 3, that IS, — 7 T ~ 
p p q q p q 


3 

2' 


Then by Lemma [TS] one has Bi{s,t,T) < 4||7’||p,||a||g(t — s). Similarly, one can establish the bound 
^ 2 ( 5 , t, T) < C(t — s), and hence B{s, t, T) < C(t — s). So (|3.8I) is proved. 

□ 

The lemmas that follow will be used in the proof of Theorem 12.41 
Lemma 3.5. Define 

pt itx _ 1 

At{x) = / e^^^ds = -, (3.9) 

Jo '•‘X 

Then for any 5 £ (0,1), there exists a constant c > 0 depending only on 5, such that 


|At(a;)| < c\tffsix), t £ [0,1], x £ M, 


(3.10) 


where 


fs{x) 


Ixl"^ 1 if |x| > 1; 
1 i/|x|<l. 


(3.11) 


Proof. In view of (|3.9L we have |At(x)| < J*\e'^^^\ds = t. So under the constraint t £ [0,1], we have 
|At(x)| < t < t^. On the other hand, from Lemma 2 from [25], with some constant C > 0, we have 
je“ - l| < C|x|^ (5 £ (0,1). So 


|At(x)| < I I < C|fx|'^|x| ^=Ct‘^|x|'^ \ 
jxj 

Combining this with (|3.1ip . we obtain (13.1011 . □ 

We quote Lemma 1 of [25] in a special case, convenient for our purposes. 

Lemma 3.6. Let 7 ^ < 1, 7 j + ^i+i > 1/2, and let 6 be such that 


0 ^^< 7^+7^+i 


7 









where i = 1,..., 4 (with 75 = 71 Then 


[ fs{yi-y2)fs{y2-y3)fs{y3-y4)fs{y4-yi)\yi\ '^"\y 2 \ '^^\y3\ '^^\y4\ < 00 , 

where fs{-) is as in 

Lemma [nm can be used to establish the following result. 

Lemma 3.7. The function 

H*{xi,X 2 ) := [ |Ai(xi •^u)At{x 2 -u)\\u\~^du (3.12) 

Jr 

is in L^(R^) for all (ai,a 2 ,/?) in the open region {(ai,a 2 ,/3) : ai, 02 ,/3 < 1, + /3 > 1/2, i = 1,2}. 

Proof. It suffices focus on the case where t G [0,1], otherwise a change of variable can reduce it to this case. 
We have by suitable change of variables and Lemma 13.51 that 

\\H(\\1^k 2) = [ \At{yi - ?/ 2 )At(j /2 - y3)At{y3 - ?/4)At(j/4 - yi)\\yi\~°''^\y2\~^\y3\~°'^\y4\~^dy 
Jr* 

<C [ fs(yi - y2)f5{y2 - y3)f&{y3 - y4)fs{y4 - yi)\yi\~‘^^\y 2 \~^\y 3 \~°‘^\y 4 \~^dy. 

Jr* 

Then apply Lemma 13.61 noting that 6 can be chosen arbitrarily small. □ 

Lemma 3.8. Define the function 

H^j^{xi,X 2 ) = Ai^t{xi,X 2 )\xiX 2 \~'^^‘^ [ \At{xiu)At{x 2 - u)\\u\~^A 2 .t{u) du, (3.13) 

Jr 


where 


Ai^t{xi,X2 ) 


L,{x,/T) L^{x 2 lT) , , L^{u/T) 

Li(l/T) Li(l/T) ’ i 2 (l/T)- 


Then for large enough T, we have Hfj,{xi,X 2 ) G 

Proof. By (12.101) and p.l4l) . for any e > 0 there exists C > 0, such that for T large enough, 

|^i,T(a;i,a;2)| < 1^(12:11" + |a;i|“")(|a;2|‘' + |a:2r'') 


and 

|^2,t(m)| < C(|m|" + lul"'^). 

Hence, with some constant C > 0, 


(3.14) 


(3.15) 

(3.16) 


\HIj.{xi,X2)\ < C |At(xi +'u)At(x2 - u)||'u| ^(|u|" + |u| ‘')du\xiX2\ “/^(Ixil'" + Ixil '')(|X2|" + |X2| "). 

Jr 

(3.17) 

Because by Lemma 15771 the function Hf in p.l2p is in L^(R^) for all (ai,a 2 ,/ 3 ) in an open region {(a,/?) : 
ai, a 2 , P < l,ai + P > 1/2, i = 1, 2}. By choosing e small enough, we infer that the right-hand side of p.l7l) 
is in L^(]R^), and the result follows. □ 















Lemma 3.9. Let Zxit) be as in i2.11\) . and let 


Zlr{t):= [ Ht^T{xi,X2) W{dxi)Widx2), ( 3 . 18 ) 


where 


Ht,T{xi,X2) = Ai^t{xi,X2)\xiX2\ 


At(a;i + u)/S.t{x2 — u)|m| ^A2 ,t{u) du 


( 3 . 19 ) 


Then Zxit) Z'rp(t), that is, the processes Zxit) and Z'rp{t) have the same finite-dimensional distributions. 

Proof. Using the spectral representation of A'(l) (see, e.g., [5], Chapter XI, Section 8): X(t) = /r \/f{x)W{dx), 

where W{-) is a complex Gaussian measure with Lebesgue control measure, and the diagram formula (see, 
e.g., [20], Chapter 5), we have 

Xiu)Xiv) - IE [X(u)X(u)] = [ e*(“"^+""^)v'/(xi)/(x2)lU(dxi)lU(dx2). 

By a stochastic Fubini Theorem (see [22]) Theorem 2.1) and Lemma [3.81 one can change the integration 
order to get (note that by (II.2p we have 'g{t) = g{x)dx): 


Zrit) = 


1 


r“+/5Li(l/T)L2(l/T) 

1 

■r“+/3Li(i/T)L2(i/r) 

1 

■r“+^Li(l/T)L2(l/T) 


'R 2 

'R 2 

'R 2 


\/f{xi)f{x 2 ) 

Vf{Xl)f{x 2 ) 


nTt pTt 


/o ^0 JR 
r pTt 


e^(^-A^g{w)dw W{dxi)W{dx2) 


pTt 


\w\-f^L{w)dw Widxi)W{dX2) 


Vf{^i)fix2) / Axtixiw)ATt{x 2 - w)\w\ ^L2{w)dw W{dxi)W{dx2). 


Now we use the change of variables w —>■ u/T, xi —?> xijT, X2 —t X2/T, where the latter two change of 
variables are subject to the rule W{dx/T) = T~^/'^W{dx) (see, e.g., [4], Proposition 4.2), to obtain 


Zrit) 


r“+/5Li(l/T)L2(l/T) 
/ ^/f{x,/T)f{x2/T) 

aR 2 


At(xi + u)At(a:2 — u)|w/T| ^L2(w/T)TdwT ^lU((ixi)lU((ix2). 

(3.20) 


Taking into account the equality fixjT) = \x/T\ °‘Li{x/T) and equations in p.l4p . we see that the right 
hand side of (I3.20p coincides with (I3.18p . This completes the proof. □ 

The lemmas that follow will be used to establish tightness in the space Cp, 1] in Theorem 12.41 

Lemma 3.10. Let 6 be a fixed number within the range (0, (a + /3)/2), and let Zxit) be as in S2.11\) . Then 
for all 0 < s < t < 1 and T large enough, there exists a constant C > 0, such that 

IE [\ZT{t) - ZTis)f] < C{t - .sf^. (3.21) 

The same estimate also holds for the corresponding limiting process Z(t) defined by S2.12\) . li2.1S\) . 

Proof. First, in view of Lemma 15^ we have IE [\ZT{t) — Zt{s)\‘^] = IE [\Zlp{t) — Zy(s)p]. Next, using the 
linearity of the multiple stochastic integral, we can write 

Z'rp{t) - Z'j.{s) = [ Hs^t,T{xi,X 2 )W{dxi)W{dx 2 ), 

JR 2 
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where 


Hs,t,T{.Xi,X2) = Ai^t{xi,X2)\xiX2\ [Aj (xi + u) A* (x2 - It) - A* (cci + m) (x2 - u)] |u| ^A2,T{u)du. 

JR. 

(3.22) 

The term in the brackets of the integrand in (13.221) can be rewritten as follows: 

At(a;i + u)At{x2 - u) - As{xi + u)As(x2 - u) 

Jo Jo Jo Jo 

PS pt pt PS pt pt 

= / dwi / dw2 ... + / dwi / dw2 ... + / dwi / dw2 ■ ■ ■ 

J 0 J s J s -/O J s J s 

=As{xi + u)At-s{x2 -u) + At-s{xi + u)As{x2 -u) + At-s{xi + u)At-s{x2 - u). 

Now we apply Lemma 13.51 to get 

|At(xi + u)At{x 2 - u)- As{xi + u)As{x 2 - w)| < C[s^{t - sY + {t- sYs^ + {t- sYYfsixi + u)fs(x 2 - u) 

<C{t-sYf5ixi+u)fsix2-u), (3.23) 

where the last inequality follows because 0 < < 1 and 0 < (t — sY < 1- 

Next, using formula (4.5') of [20], (13.221) and (13.231) . we can write 

]E[\ZT(t) - ZTis)\Y = ||i?s,t.r|li2™2', < Clt - [ dxidX2Ai^Tixi,X2Y\xiX2\~°'x 

/ dUidU2fs{xi + Ui)fs{x2 - Ui)fs{-Xi + U2)fs{-X2 - U2 )\ui\~^\u2\~'^A2^t{ui)A2^t{u2) 
dR2 

<(71^ — sp"* / dyidy2dyody4:Ai^T{yi,yd,YA2^T{y2)A2^T{yi)^ 

JRO 

fs{yi - y2)f5iy2 - yo)fs{yo - yYfsiy^ - 2/i)|yir“l2/2r'^|j/3r“|y4r'^, (3.24) 

where we have applied the change of variables: yi = xi, y 2 = —ui, yo = —X 2 , j /4 = U 2 . 

Since by assumption a < 1, j3 < 1 and a + /3 > 1/2, and the exponent e in ()3.15|) and (13.161) can be 
chosen arbitrarily small, for a fixed 6 satisfying 0 < 6 < (a + /3)/2, we can apply Lemma 1 of [5^ to conclude 
that the integral 

/ dyAi^Tiyi,y3YA2,T{y2)A2,T{y4)fs{yi - y2)fs{y2 - y3)fs{y3 - yYfsiyi - yi)\yi\~°‘\y2\~^\y3\~°‘\y4\~^ 


is bounded for sufficiently large T, which in view of (13.241) implies (13.211) . The proof for ZtY) is thus 
complete. The proof for Z{t) is similar and so we omit the details. □ 


4 Proof of Main Results 

Proof of Theorem \2.1[ By Lemma [321 for any 0 < ti < ... < and constants ci,..., c„, we have 


lim Var 

T —>-oo 


n 

(QT{tj) - Lritj)^ 

1=1 


= 0 . 


Therefore the convergence of finite-dimensional distributions of QtY) to that of Brownian motion aB(t) 
follows from Lemma [3.11 with /y(-) given in (13.71) and the Cramer-Wold Device. □ 


10 


















Proof of Theorem \2.‘A In view of the well-known Prokhorov’s Theorem (see, e.g., [3], p. 58), to prove the 
theorem, we need to show convergence of finite-dimensional distributions and tightness. The former has 
been established in Theorem l2.ll To prove tightness, observe that by Lemma [3131 and the hypercontractivity 
inequality of the multiple Wiener-Ito integrals (see [3D], Corollary 5.6), for any T > 0 and 0 < s < t < 1, 
there exists a constant C > 0 to satisfy 


IE 


IQrit) - Qt(s)|^] < C 2 (lE [iQT(t) - Qt(s)|"])" < C(t - s) 


(4.1) 


Now the tightness of the family of measures generated by the processes {Qrit) : T > 0} in Cp, 1] follows 
from Lemma 5.1 of m- □ 


Proof of Theorem \2.A The convergence of finite-dimensional distributions follows from Theorem 12.11 In 
fact, the assumptions on / and g in Theorem 12.31 imnlv the conditions (|2.2I) and (|2.3I) in Theorem 12.11 fsee 
the proof of Theorem 5 in nni). The tightness can be shown similarly as in the proof of Theorem 1 2.2 1 □ 


Proof of Theorem 12-41 in the proof of Theorem 12.21 we need to show convergence of finite-dimensional 
distributions and tightness. We first prove the convergence of finite-dimensional distributions, that is, 
Zrit) Z(t) as T —>• oo, where Zrit) and Z(t) are defined by (|2.11l) and (12.121) . respectively. 

By Lemma [3.91 the process Zrit) defined in (12.111) has the same finite-dimensional distributions as the 
process Zf,{t) defined in (I3.18p . Therefore, taking into account the linearity of multiple Wiener-Ito integral, 

and applying Cramer-Wold device, to prove Zxit) Z{t), it is enough to show that as T —)> oo, 

Ht,T{xi,X 2 ) ^ Ht{xi,X 2 ) in L^(R^), (4.2) 

where Ht{xi,X 2 ) and Ht^T{xi,X 2 ) are as in (12.131) and (I3.19L respectively. 

First, we show pointwise convergence for a.e. (xi,X 2 ) S that is. 


Ht,Tixi,X 2 ) = Ai^Tixi,X 2 )\xiX 2 \ / At{xi + u)At(x 2 - u)\u\ ^A 2 ^T{u)du 

Jr 

Ht{xi,X 2 ) = \xiX 2 \~°‘^‘^ j At{xi + u)At{x 2 — u)\u\~^du as T —>• oo. 
JR 


(4.3) 

(4.4) 


Because Li{x) is a slowly varying function, we have Ai^t(xi,X 2 ) —>■ 1 as T —?> oo, where Ai^t is as in (13.141) . 
To show that the integral in (14.31) converges to the integral in (14.4L note first that by (13.141) . A 2 ^t{u) —>■ 1 as 
T —>■ 00 because L 2 {x) is a slowly varying function. Hence one only needs to bound the integrand properly 
and apply the Dominated Convergence Theorem. To this end, observe that by p.l 6 l) for T large enough, we 
have 


gT{u]Xi,X 2 ) ■ = |At(a:i -I- u)||At(a ::2 - u)\\u\~^A 2 ,t{u) (4.5) 

< C'|At(a;i -f u)||At(a:2 - u)||u|“'^(|'u|'' -f lul"") := ^^(m;X i,X 2 ). (4.6) 

By choosing e small enough, using Fubini Theorem and Lemma 13.71 we conclude that ge{- ;xi,a; 2 ) G L^(K.) 
for a.e. (xi,a; 2 ) £ R^. Now (14.21) follows from (13.171) and the Dominated Convergence Theorem. 

To prove tightness, first observe that by the hypercontractivity inequality of the multiple Wiener-Ito 
integrals (see [3D], Corollary 5.6) and Lemma [3.101 for T large enough and for any 0 < s < t < 1, there 
exists a constant C > 0 to satisfy 

W.[\ZT{t) - Zt{s)\^] <C 2 {m[\ZT{t)-ZT{s)\^]f <C\t-s\^\ (4.7) 

where d is a fixed number within the range 0 < 45 < 2(a -|- jd). Since by assumption a + fd > 1/2, we can 
choose 5 to satisfy 45 > 1. Inequalities similar to (14.71) hold also for the limit process Z{t). 

In view of (14.71) and a similar inequality for Z{t), it follows from Kolmogorov’s criterion (see, e.g., [3] 
Theorem 8.1(1)) that the processes Z'r(t) and Z{t) admit continuous versions when T is large enough. 

Now the tightness of the family of measures generated by the processes {^^(t) : T > 0} in Cp, 1] follows 
from Lemma 5.1 of m- Theorem 12.41 is proved. □ 
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